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ABSTRACT 


Introduction: Developing and optimizing isocratic HPLC methods is a difficult procedure that requires simultaneous determina- 
tion of several factors. 


Objective: To develop a simple, accurate, rapid, and precise chemometric-assisted RP-HPLC method for the simultaneous 
estimation of Paracetamol (PARA) and Ibuprofen (IBU) in bulk and in pharmaceutical formulation. 


Methods: The method involves the design of experiment (DOE) approach for the optimization of the mobile phase by keeping in 
concern about the dependent variables (flow rate, methanol, and pH) and their consequence was perceived over the retention 
time (Rt) and resolution of drugs. 

Results: Based on the trial runs, the optimized condition generated by the software was used for the separation of PARA and 
IBU on Thermo C18 Column (150 mm x 4.6 mm i.d., 5 um). The chromatographic separation was achieved by using the op- 
timized mobile phase consisting of methanol and sodium dihydrogen phosphate buffer (88.3: 11.7 % v/v), adjusted to pH 2.5, 
delivered at a flow rate of 1.2 mL/min and the detection was carried out at 222 nm. The retention time (Rt) was found to be 2.241 
minutes for PARA and 3.535 minutes for IBU. The developed method was validated as per ICH guidelines. The calibration curves 
of PARA and IBU were found to be linear in the range of 6.5-39 ug/mL and 8-48 ug/mL, respectively. The recoveries for the two 
compounds were close to 100% at various levels. 

Conclusion: This method can be successfully applied for routine quality control analysis of these drugs in bulk and pharmaceuti- 
cal formulations. 
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INTRODUCTION is a well-known NSAID preferably act by the inhibition 


of COX-2.^ 
Paracetamol (PARA), chemically known as N-(4-hydroxy- 


phenyl) acetamide (Figure 1A) is the widely employed 
non-steroidal anti-inflammatory drug (NSAID) that acts 
by inhibition of cyclooxygenase-1/2 (COX-1/COX-2) en- 
zyme which prevents the metabolism of arachidonic acid 
to prostaglandins H, and E, thereby lower the hypothalam- 


While looking at some imperative medical databases such 
as Google Scholar, Pubmed, etc. for quality analysis (esti- 
mation) of these two drugs (PARA and IBP) combination 
through detailed literature survey, it was observed that cer- 
tain methods like Ultraviolet-Visible (UV-Vis),>"'° High-Per- 


ic set-point to reduced fever.'? Ibuprofen (IBP), an aryl 
acetic acid derivative that is chemically known as (RS)- 
2-(4-(2-methylpropyl)phenyl)propanoic acid (Figure 1B) 
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was noticed that no article was available in the literature re- 
garding the “chemometric approach (design of experiment)” 
used for the simultaneous estimation of PARA and IBP using 
RP-HPLC. As the chemometric methods are well-known to 
rationally offer the best logical outcomes, hence, this ap- 
proach needs to be developed. 





Developing and optimizing isocratic HPLC methods 1s a dif- 
ficult procedure that requires simultaneous determination 
of several factors, to optimize more than one response at a 
time and to gain information on the sensitivity of chroma- 
tographic factors along with their interaction effects on the 
separation characteristics. Hence, the chemometric method 
which includes the robust applications of factorial design 
and response surface methodology.???! 





This research investigation describes the development and 
validation of a fast, easy, and sturdy RP-HPLC method for 
the simultaneous estimation of PARA and IBP using chemo- 
metric-design of experiment (DOE) approach. 


EXPERIMENTAL 


Chemicals 

The working standards of PARA and IBP were received as 
generous gift samples from Himalaya Meditek Pvt. Ltd., 
Dehradun, Uttarakhand. Disodium hydrogen phosphate (AR 
grade) and methanol (HPLC grade) were procured from 
Loba Chemie Ltd., Mumbai, Maharashtra. Other analytical 
grade chemicals, reagents, and solvents were obtained from 
HiMedia Pvt. Ltd., Mumbai, Maharashtra. HPLC grade wa- 
ter was purchased from the ELGA system. Combiflam® tab- 
lets, manufactured by Sanofi India Ltd., Gujarat containing 
PARA I.P. 325 mg and IBU I.P. 400 mg were acquired from 
a local pharmacy shop. 





Instrumentation 

A JASCO (Japan) HPLC system, equipped with PU 2080 
Plus solvent delivery system (pump), MD 2010 plus PDA 
(detector), and Rheodyne injector with 20 uL loop volume 
was employed. The data were obtained and processed by uti- 
lizing the Borwin software v.1.5. 


Chromatographic conditions 

The chromatographic separation was achieved utilizing the 
Thermo? C,, column (150 mm x 4.6 mm, i.d., 5 um). The 
composition of the mobile phase consisted of methanol and 
sodium dihydrogen phosphate buffer (88.3: 11.7 9o v/v; pH 
2.5 adjusted with 10% orthophosphoric acid), delivered at a 
flow rate of 1.2 mL/min. The column effluents were moni- 
tored at 222 nm. 


Preparation of phosphate buffer solution 
The phosphate buffer solution was prepared by dissolving 90 
mg of disodium hydrogen phosphate in 100 mL of distilled 


water to produce 10 uM concentration. The pH was adjusted 
to 2.5 with 10% orthophosphoric acid. 


Preparation of standard stock solution 

The stock solution of PARA and IBP were prepared by dis- 
solving accurately weighed 10 mg of standard drug in 10 mL 
of methanol, separately to produce 1000 ug/mL concentra- 
tion each. Further, 1 mL of this solution was diluted to 10 mL 
to produce 100 ug/mL concentration each. 





Preparation of working stock solution 

The working solutions were prepared from the standard stock 
solution (100 ug/mL) by appropriate dilution in the mobile 
phase to obtain the final concentration of 6.5-39 ug/mL for 
PARA and 8-48 ug/mL for IBP, respectively. 





Optimization 

The optimization of the mobile phase conditions was per- 
formed as per the experimental design by employing a three 
factor-two level (3?) response surface method of Central 
Composite Design (CCD) using the Design-Expert software 
v.9.0.6 (Stat-Ease Inc., Minneapolis, USA). Methanol vol- 
ume (A), flow rate (B), and pH (C) were selected as the in- 
dependent variables, while the retention time of IBP (tR2) 
and resolution (Rs2) were taken as responses. The response 
surface analyses were performed to recognize the effect of 
different independent variables on the observed responses 
(Table 1). A total of 20 runs were generated by the Design- 
Expert software and all the experiments were carried out in 
a randomized order to reduce the responses of uncontrolled 
variables that may bring in unfairness on the measurements. 
The responses were thoroughly evaluated statistically em- 
ploying the ANOVA process and were observed to be sta- 
tistically significant. Additionally, the most favourable state 
was chosen through the numerical optimization procedure 
employing the desirability function. As the P-value was 
found to be «0.05, the quadratic model was used. The quad- 
ratic mathematical model for independent factors 1s depicted 
by the equation below: 























Y = By + BX, t B,X, + B,X, + B,;X,X,  B,,X,X4 + B23X2X3 + By Xi + B2;X2 
+ Bat 


Global Optimization 

In the present study, to optimize the two responses with dif- 
ferent targets, Derringer’s desirability function was used. 
The Derringer’s desirability function (D) is defined as the 
geometric mean, weight, or otherwise, of the individual de- 
sirability functions. The expression that defines the Derrin- 
ger’s desirability function is given by the equation below. 


D = [d?*x d? xa? x... xa" |" 


Where, n- the number of responses and Pn- the weight of 
the responses 
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The weight of the response is the relative importance of each 
of the individual functions di. The relative importance pi is 
a comparative scale for weighting each of the resulting di in 
the overall desirability product and it varies from the least 
important (pi = 0.1) to the most important (pi = 10). The 
desirability function (D) can take values from 0 to 1. The 
weights can range from 0.1 to 10. The weights <1 give less 
significance to the desired objective whereas the weights >1 
described better significance to the desired objective. In the 
current research, pi values were set at 1 for all the three re- 
sponses. A value of D close to 1 indicates that the combina- 
tion of the different criteria is matched in a global optimum. 
The percentage of prediction error was calculated by the 
equation below. 








Experimental value — Predicted value 


Prediction Error = x 100 


Predicted value 

Method Validation 

The developed RP-HPLC method was properly corroborated 
in agreement with the Q2A and Q2B guidelines of the Inter- 
national Council for Harmonization of Technical Require- 
ments for Pharmaceuticals for Human Use (ICH), in assent 
with the United States Pharmacopoeia (USP), and with the 
guidance of the United States Food and Drugs Administra- 
tion (USFDA). 





Specificity 

The degree of specificity was determined by injecting the 
blank solution and placebo (synthetic mixture) solution, and 
further comparing the peaks in the sample solution.? 


Linearity 

For the evaluation of the system linearity, five concentrations 
(80%, 90%, 100%, 110%, and 120%) of the target analytical 
concentration were taken. The solutions were carefully pre- 
pared in the presence of a diluent and the equivalent volumes 
were injected under the particular chromatographic condi- 
tion. The linearity equation was derived based on the aver- 
age areas versus the concentration volume and the regression 
coefficient value (r?) was reported.? 





Accuracy 

The accuracy (recovery) of this RP-HPLC method was esti- 
mated by spiking the standard drug in the concentrations of 
50%, 100%, and 150% of the target concentration (standard 
addition method). The process was accomplished in a tripli- 
cate way and the acquired mean data were articulated in the 
form of % recovery + confidence interval with calculated % 
relative error, based on the specific concentrations.” 


Precision 
The precision (variability) of this RP-HPLC method was de- 
termined by spiking the standard drug in the concentrations 


of 80%, 100%, and 120% of the target concentration (stand- 
ard addition method), three-times in a single day (intra-day 
variability) and three-times on three different days (inter-day 
variability). The level of precision was determined through 
the acquired relative standard deviation (RSD) data.” 


Robustness 

The deliberate alteration in the systems suitability parame- 
ters like mobile phase composition by +2% v/v (1.e., 90.3:9.7 
% v/v and 86.3:13.7 % v/v) and detection wavelength by 
+2°C (1.e., 224 nm and 220 nm) were methodically studied 
keeping other factors constant.”° 





System suitability 

The reproducibility characteristics of the system were esti- 
mated by spiking the standard solution five-times and deter- 
mining the parameters like tailing factor (TF), retention time 
(Rt), theoretical plates (TP), and peak area.” 


Limit of detection (LOD) 

Limit of detection (LOD) may be defined as the lowest con- 
centration that can be detected by the method but not always 
necessary to quantify in the exact value. The LOD was deter- 
mined by the following formula:?? 


LOD - 3.3 (o/S) 


where o = standard deviation of response and S = slope of 
the calibration curve. The slope S may be estimated from the 
calibration curve of the analyte. 


Limit of quantification (LOQ) 

Limit of quantification (LOQ) may be defined as the low- 
est concentration that can be quantified dependably with a 
particular level of accuracy and precision. The LOQ is deter- 
mined by the following formula:” 


LOQ = 10 (o/S) 


where o = standard deviation of response and S = slope of 
the calibration curve. The slope S may be estimated from the 
calibration curve of the analyte. 


Assay of Combiflam® Tablets 


20 tablets of Combiflam? were accurately weighed and fine- 
ly powdered. The powder equivalent to 10 mg of PARA and 
12.30 mg of IBP was taken and transferred to a 10 mL volu- 
metric flask. The content was diluted to 10 mL with metha- 
nol and the solution was sonicated for 10 min. The solution 
was then filtered through Whatman filter paper no. 41. 1 mL 
of the filtrate solution was diluted to 10 mL with the mobile 
phase. Further 1.3 mL of this solution was diluted to 10 mL 
with methanol to achieve a final concentration of 13 ug/mL 
and 16 ug/mL, respectively of PARA and IBP. This proce- 
dure was repeated 6 times for the tablet formulation. 


Int J Cur Res Rev | Vol 12 + Issue 23 - December 2020 100 


Kanthale et al.: Development of a rapid chemometric reverse phase hplc method for the simultaneous estimation of paracetamol and ibuprofen 


Forced degradation studies 


Alkali degradation studies 

0.5 mg equivalent of the drug was taken in a 100 mL volu- 
metric flask and diluent (methanol) was added to exactly the 
half, sonication was done for 15 min duration, and further, 
the volume was made up to 100 mL. The content was stirred 
using the magnetic stirrer for 30 min duration and then at 
3,000 rpm centrifugation was done for 5 min. 5 mL of the 
above content was taken in a volumetric flask, 10 mL of 0.1 
N NaOH solution was added, and the above content was 
heated at 80°C temperature for the duration of 8 hrs under 
dark condition. The content was neutralized with the same 
quantity of HCI and the volume was made up to 100 mL with 
the mobile phase. The final content was filtered employing 
the 0.45 um pore size nylon membrane.?? 


Acid degradation studies 

0.5 mg equivalent of the drug was taken in a 100 mL volu- 
metric flask and diluent (methanol) was added to exactly the 
half, sonication was done for 15 min duration, and further, 
the volume was made up to 100 mL. The content was stirred 
using the magnetic stirrer for 30 min duration and then at 
3,000 rpm centrifugation was done for 5 min. 5 mL of the 
above content was taken in a volumetric flask, 10 mL of 0.1 
N HCI solution was added, and the above content was heated 
at 80°C temperature for the duration of 8 hrs under dark con- 
dition. The content was neutralized with the same quantity 
of NaOH and the volume was made up to 100 mL with the 
mobile phase. The final content was filtered employing the 
0.45 um pore size nylon membrane?! 


Oxidation degradation studies 

The oxidative stress (hydrogen peroxide method) tendered 
to the drug was determined by dissolving 0.5 mg equivalent 
weight in 10 mL of H,O, (3% v/v) in a volumetric flask and 
further boiling for the duration of 1 hr. The above content 
was kept at the room temperature to commence the process 
of degradation. Then, the content was suitably diluted, fur- 
ther sonicated, and the volume was carefully made up to 100 
mL volume. The final content was filtered employing the 
0.45 um pore size nylon membrane.” 


Photostability studies 

In photostability studies, the compound absorbs the photon 
and leads to the overlapping of the absorption band which 
results in valence electron reaching the excited state. 0.5 mg 
equivalent of the drug was taken in a petri dish and exposed 
to sunlight for the duration of 12 hrs for a single day. The ex- 
posed content was further transferred in a 100 mL volumetric 
flask, the mobile phase was added, sonication was performed 
for 15 min, and the volume was made up to 100 mL. The 
final content was filtered employing the 0.45 um pore size 
nylon membrane. The radiations of 254 nm wavelength com- 


menced the degradation process and the primary degradation 
products were determined by RP-HPLC.” 


Thermal studies 

In the thermolytic study, the drugs instigate degradation via 
assorted mechanisms under elevated temperatures (>70°C) 
in the solution states. The study was performed by taking 0.5 
mg content equivalent to the drug, exposed to heat at 8041°C 
for 1 hr, and evaluated for stability under the stress environ- 
ment. 10 mL of methanol was added, sonicated for 15 min, 
and volume was made up to 100 mL with the mobile phase. 
The sample solution was then stirred for 30 min employing 
the magnetic stirrer, and filtered through a 0.45 um pore size 
nylon membrane. 20 uL of drug solution was injected into 
HPLC. 








RESULTS AND DISCUSSION 


Optimization 

DoE is an approach that enables researchers to assess the 
consequences of potential interactions of several variables 
on productivity concurrently using a restricted number of 
experiments. From the data generated in preliminary studies 
while developing the separation for PARA and IBU, limits 
for experimental levels were eventually recognized. 


The statistical parameters obtained from ANOVA for the re- 
duced model are given in Table 2. To get a more realistic 
model, unimportant terms with corresponding p-value 70.05 
were removed during the backward elimination process. The 
adjusted R? values were well within the satisfactory limits 
of R? > 0.80, publicized that the experimental data are in 
good fit with the second-order polynomial equations. Since 
P-value is « 0.05, reveals that all the reduced models are sig- 
nificant. In this study, the signal (response) to noise (devia- 
tion) ratio was found to be in the range of 9.04-17.44 (ratio 
greater than 4 1s desirable), suggestive of an adequate signal 
to noise ratio and therefore the model is significant for the 
separation process. The %CV (coefficient of variation) of all 
the models were found to less than 10% revealed that all the 
models were reproducible, (model can be considered reason- 
ably reproducible if % CV is less than 10%). 








As can be seen from the regression equation (Table 2), the 
positive interaction of factors A, B, and C are statistically 
important for tR, and Rs,. Altering the factor A and B, low to 
high level strongly effect (decreasing order) the tR, and Rs,. 
Then altering factor C show a slight influence tR and Rs, 
(lower level C and A show highest resolution). Therefore, 
this study discloses that intensifying factors A, B, and C will 
reduce tR, and Rs,. Hence, a high level of factors A, B, and 
C will give a shorter run time. 


Perturbation plots and 3D response surface plots are depicted 
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for the predicted models, to give an improved understanding 
of the investigated method. This type of plot represents the 
effect of an independent factor on a specific response with all 
other factors assumed constant at a reference point. A steep- 
est slope or curvature represents the sensitiveness of the re- 
sponse to specific factors. Figure 2 displays the percentage 
of methanol in the mobile phase (A) is the most important 
effect on tR, and Rs, followed by factor B for tR, and factor 
C for Rs, The factor A and B are important for decreasing the 
tR, and Rs, but factor C is also having a little effect on Rs,. 
So these three factors (A, B, and C) were important for the 
shorter run time. 


The criteria for the optimization of each response are shown 
in Table 3. As can be seen under the criteria for optimization, 
the responses tR, and Rs, were minimized and targeted to 
shorten the analysis time and give better resolution respec- 
tively. Following the conditions and restrictions above, the 
optimization procedure was carried out. The response sur- 
face attained for the global desirability function 1s presented 
in Figure 3. It can be concluded that there was a set of co- 
ordinates producing high desirability value (D= 0.862) were 
methanol concentration (A) of 88.30 %, pH 2.5 and flow rate 
of 1.20 mL/min. The prediction competence of the model 
was substantiated by executing the experiment under the 
most favourable state and the respective chromatogram. To 
attain the obviousness of the suggested model, the assertion 
between the experimental responses and predicted responses 
for the envisaged criteria is shown in Table 4. 











Chromatographic condition 

The chromatographic separation was achieved using Ther- 
mo^ C, column (150 mm x 4.6 mm, i.d., 5 um). The opti- 
mum composition of the mobile phase consisted of methanol 
and sodium dihydrogen phosphate buffer (88.3: 11.7 96 v/v; 
pH 2.5 adjusted with 10% orthophosphoric acid) at a flow 
rate of 1.2 mL/min, which resolved the degraded peak from 
the standard peak with Rt values of 2.241 min for PARA and 
3.535 min for IBU (Figure 4). The UV spectra showed that 
the drug absorbed appreciably at 222 nm, and therefore, this 
wavelength was selected as the detection wavelength. The 
developed method was also found to be specific, since, it 
was able to separate other excipients present in the dosage 
form of the drug. 


METHOD VALIDATION 


Specificity 

The specificity of this method was estimated where no inter- 
ference was observed. The observed peaks in the sample so- 
lution were the actual peak of PARA and IBP which showed 
that the method was specific. 


Linearity and range 

Over the desired range of 80% to 120% of the target ana- 
lytical concentration, an excellent level of linearity was ob- 
served for both the drugs. The linear regression equation for 
PARA was observed to be y = 52185x - 40419 with a re- 
gression coefficient value (r?) of 0.996 and the linear regres- 
sion equation for PARA was perceived to be y = 120898x 
+ 727251 with a regression coefficient value (r?) of 0.995, 
which signifies a highly acceptable degree of linearity. 











Accuracy 

The accuracy characteristics of the developed method were 
computed based on the recovery data by employing the 
calibration curve where the Y-intercept and the slope play a 
critical role in estimating the 9o recovery. The result of the 
recovery study for PARA was observed to be 101.11% w/w 
(at 50% concentration), 99.81% (at 100% concentration), 
and 101.93% (at 150% concentration). The % RSD values 
were observed to be 1.05%, 1.82%, and 0.95%, respectively 
at three different concentrations; 50%, 100%, and 150% of 
the target concentration. The result of the recovery study for 
IBU was observed to be 99.81% w/w (at 50% concentra- 
tion), 101.66% (at 100% concentration), and 101.23% (at 
150% concentration). The % RSD values were observed to 
be 0.78%, 1.30%, and 1.37%, respectively at three different 
concentrations; 50%, 100%, and 150% of the target concen- 
tration. The % RSD values were found to be within the phar- 
macopoeia acceptance limit of +2% which indicated good 
accuracy of the developed analytical method. 








Precision 

Over the 50%, 75%, and 100% of the target concentration, 
the % RSD values were observed to be 1.37%, 0.98%, and 
1.54% for PARA and 1.21%, 1.40%, and 0.95% for IBU, 
respectively under intra-day variability whereas 1.48%, 
1.05%, and 1.16% for PARA and 0.84%, 1.72%, and 1.92% 
for IBU, respectively under inter-day variability which rep- 
resents that the developed chromatographic method has high 
precision value with reduced variability in the estimation of 
PARA and IBU. The % RSD values were found to be within 
the pharmacopoeia acceptance limit of +2%. 





Robustness 

After intentional alterations in the chromatographic condi- 
tions such as mobile phase composition by +2% v/v (1.e., 
90.3:9.7 % v/v and 86.3:13.7 % v/v) and detection wave- 
length by +2°C (1.e., 224 nm and 220 nm), minute changes 
in the retention time (20.2) were observed. The % RSD val- 
ues were observed to be within the pharmacopoeia accept- 
ance limit of +2% which concluded that the method is robust 
enough to detect PARA and IBU even after modifications in 
the chromatographic environment. 
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Systems suitability parameters 

The system suitability parameters revealed that the devel- 
oped analytical method has reproducible attributes, capa- 
bility for routine analysis, and competent enough to be ap- 
plied. The average retention time (Rt) of PARA was noticed 
at 2.227 min and IBU was perceived at 3.533 min (Table 
5). The number of theoretical plates was recognized to be 
2303.029 for PARA and 2278.903 for IBU which are above 
the prescribed USP limit of 2,000 that may be translated into 
higher column efficacy, better resolution, and superior sepa- 
ration of the method. The high value of peak area may be 
correlated with the theoretical separation ability of the de- 
veloped method. The tailing factor was perceived to be 1.62 
for PARA and 1.41 for IBU which represented a good sym- 
metry where the degree of asymmetric factor is also equal 
to 1, an ideal form of Gaussian peak where both are factors 
hold equal magnitude. The method has been found to fulfil 
the minimum requirements as prescribed in the monographs 
of USP and has perspectives to express reproducible results. 














Forced degradation studies 

After alkaline treatment, PARA and IBU showed no peak 
of degradation product (Figure 5A). The percent recovery 
of PARA and IBU was 86.45% and 91.57%, respectively. 
After acid hydrolysis, PARA showed one peak of degraded 
product at Rt value of 2.824 min with 65.81 % recovery. On 
the other hand, IBU was 80.42 % recovered, with no peak of 
degradant (Figure 5B). In the oxidative condition, the per- 
cent recovery obtained for PARA was 70.46 % and IBU was 
96.64%, with no peak of degradant (Figure 5C). In the dry 
heat degradation condition, the percent recovery obtained for 
PARA was 95.83% and IBU was 92.10%, with no peak of 
degradant (Figure 5D). After the photo-degradation study 
for UV light and Fluorescence light, PARA was 98.80 % and 
IBU was 98.58 % recovered, with no peak of degradant (Fig- 
ure 5E). 





Analysis of Marketed Formulation 

The method was successfully used for the determination 
of PARA and IBU in pharmaceutical formulation Combi- 
flam®. The assay results for these NSAID tablets using the 
proposed method showed a recovery of 100.69% for PARA 
and 100.82% for IBU which indicated the suitability of the 
method to determine both the contents in the tablet formula- 
tions, without detecting the excipients (Table 6). It can be 
concluded that the method separates the drug from their deg- 
radation products and it may be employed for the routine 
analysis of PARA and IBU in tablets. 


CONCLUSION 


This methodology represents a well-organized and easy 
complete approach to conclude the problem of probing op- 


timum RP-HPLC conditions. This investigation also proved 
that the chromatographic techniques coupled with chemo- 
metric tools, afford useful information regarding separation 
and elution time and hence make this combined technique a 
powerful analytical tool. The validation study also proved 
that the developed and optimized method was specific, suit- 
able, linear, accurate, precise, and robust for the estimation 
of PARA and IBP. Therefore, this HPLC method can be used 
as a routine quality control analysis in a pharmaceutical en- 
vironment. 
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Table 1: Central composite design (CCD) consisting of experiments for the study of three factors with two 


responses. 
FACTOR 1 FACTOR 2 
(A): METHANOL (B): FLOW RATE 
(%) (ml/min) 
1 81.60 1.00 
2 85.00 0.80 
3 85.00 1.20 
4 85.00 0.80 
5 85.00 1.20 
6 90.00 1.00 
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FACTOR 3 RESPONSE 1 RESPONSE 2 
(C): pH Retention time Resolution 

(tR.) (Rs,) 
3.00 8.640 4.220 
2.50 9.603 5.579 
2.50 6.183 3.109 
3.50 5.660 1.222 
3.50 4-927 2.335 
3.00 4.220 1.240 
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Table 1: (Continued) 





FACTOR 1 FACTOR 2 FACTOR 3 RESPONSE 1 RESPONSE 2 
(A): METHANOL (B): FLOW RATE (C): pH Retention time Resolution 
(%) (ml/min) (tR) (Rs) 
7 90.00 1.00 3.84 4.626 1.757 
8 90.00 1.00 3.00 4.273 1.310 
9 90.00 1.00 3.00 4.210 1.242 
10 90.00 1.00 3.00 4.223 1.274 
11 90.00 1.33 3.00 3177 1.848 
12 90.00 1.00 3.00 4.270 1.250 
13 90.00 1.00 3.00 4.225 1.244 
14 90.00 0.66 3.00 6.397 1.974 
15 90.00 1.00 2.15 4-927 1.801 
16 95.00 0.80 3.50 4.790 0.838 
17 95.00 1.20 2.50 3.593 1.811 
18 95.00 1.20 3.50 3.255 0.748 
19 95.00 0.80 2.50 5.333 1.613 
20 98.40 1.00 3.00 3.541 0.870 


Table 2: Reduced response model and statistical parameters obtained from ANOVA (After backward elimina- 
tion). 


RESPONSE tR, Rs, 


Regression Equation +4.23 - 1.29À - 0.91B - 0.45C + 0.053AB + 0.48AC 147 - 0.92A - 0.13B - 0.49C + 0.14AB + 0.45AC + 
+ 0.31BC + 0.68A? + 0.22B? + 0.21C? 0.45BC + 0.54A? + 0.24B?+ 0.20C’ 

Adjusted R? 0.909 0.856 

Model p-value < 0.0001 «0.0001 

Adequate Precision 17.447 9.049 

% CV 9.73 8.46 


Table 3: Criteria for optimization of individual responses. 





RESPONSE LOWER LIMIT UPPER LIMIT CRITERIA-I 

Goal Importance 
tR, 3177 9.180 Minimize 3 
Rs, 0.748 5.579 Target = 2 3 


Table 4: The comparison of observed and predicted values under optimum condition. 


OPTIMUM CONDITIONS METHANOL (%) FLOW RATE (mL/min) pH tR, Rs, 
Predictive 88.30 12 2.5 3.547 2.000 
Experimental 88.30 1.2 2.5 3.533 2.063 
Percent Error (96) - - - 0.396 3.15 


Desirability value (D) = 0.862 
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Table 5: Validation results of developed and optimized reverse phase chromatographic method. 





PARAMETERS PARACETAMOL IBUPROFEN 
Specificity No interference No interference 
System suitability 
a) RT 2:227 min 3.533 min 
b) Theoretical plates 2303.029 2278.903 
C) Asymmetric factor 1.62 1.41 
d) Resolution - 2.063 
Range 6.5- 39 pg/mL 8-48 pg/mL 
Linearity y = 52185x - 40419 y = 120898x + 727251 
R? = 0.996 R? = 0.995 
Accuracy % Recovery % RSD % Recovery % RSD 
a) 5o % 101.11 1.05 99.81 0.78 
b) 100 96 99.81 1.82 101.66 1.30 
c) 150 % 101.93 0.95 101.23 1.57 
Precision- Intra day Concentration % RSD Concentration % RSD 
a) 13 L37 a) 16 1.21 
b) 19.5 0.98 b) 24 1.40 
c) 26 1.54 c) 32 0.95 
Precision- Inter day Concentration % RSD Concentration % RSD 
a) 13 1.48 a) 16 0.84 
b) 19.5 1.05 b) 24 1.72 
C) 26 1.16 c) 32 1.92 
Robustness Detection value % RSD Detection value % RSD 
a) Wavelength 
UEM 224 0.74 224 1.79 
ici 220 1.25 220 1.56 
b) Mobile phase 
(0) 
Y M 90.57 1.83 90.57 1.48 
“2 87.03 1.49 87.03 1.72 
LOD 2.06 1.56 
LOQ 6.26 4.75 


Table 6: Assay result for PARA and IBP in Combiflam tablet by proposed methods. 





PARACETAMOL IBUPROFEN 
Actual Peakarea AmountRe- %Recovery Actual Peakarea Amount Recov- % Recovery 
amount covered (ug/ amount ered (pg/mL) 
ml) 

1 13 654834.9 13:32 102.48 16 2697852.0 16.29 101.87 
2 13 640342.6 13.04 100.34 16 2689896.9 16.23 101.46 
2 13 644723.1 13.12 100.99 16 2699887.6 16.31 101.97 
4 13 640045.8 13.03 100.30 16 2699999.0 16.31 101.98 
5 13 639987.4 13.03 100.29 16 2629896.8 15.73 98.36 
6 13 636342.8 12.96 99.75 16 2647690.9 15.88 99.27 
Mean 100.69 100.82 
SD 0.959 1.590 
jmeb 0.952 1.577 
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Figure 1: Structure of (A) Paracetamol (B) Ibuprofen. 
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Figure 2: 3D response surface plots: (A) methanol, (B) flow rate, (C) pH, and (D) perturbation plots. 


Desirability 





A:MeOH 1 
3:flow rate 1 
C:pH 1 
t2 0.771767 
res 1 
Combined 0.862396 
0.000 0.250 0.500 0.750 1.000 


Figure 3: Depiction of global desirability function. 
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Figure 5: Chromatograms of forced degradation investigation under conditions of: (A) Alkaline treatment, (B) Acid treatment, (C) 
Oxidative stress, (D) Dry-heat environment, and (E) Photo-degradation. 
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